Mon. Not. R. Astron. Soc. 000. ITHT91 (2009) Printed 4 January 2010 (MN I^TeX style file v2.2) 



What determines satellite galaxy disruption? 



Andrew R. Wetzel 1 , Martin White 1 ' 2 

1 Department of Astronomy, University of California, Berkeley, CA 94720, USA 
2 Department of Physics, University of California, Berkeley, CA 94720, USA 



o 

(N 



O 

U 

6 



> 

(N 
O 

q 

o 

o\ 
o 



>< 



July 2009 



ABSTRACT 

In hierarchical structure formation, dark matter halos that merge with larger halos can 
persist as subhalos. These subhalos are likely hosts of visible galaxies. While the dense 
halo environment rapidly strips subhalos of their dark mass, the compact luminous 
material can remain intact for some time, making the correspondence of galaxies with 
severely stripped subhalos unclear. Many galaxy evolution models assume that satellite 
galaxies eventually merge with their central galaxy, but this ignores the possibility of 
satellite tidal disruption. We use a high-resolution iV-body simulation of cosmological 
volume to explore satellite galaxy merging/disruption criteria based on dark matter 
subhalo dynamics. We explore the impact that satellite merging/disruption has on the 
Halo Occupation Distribution and radial profile of the remnants. Using abundance 
matching to assign stellar mass/luminosity to subhalos, we compare with observed 
galaxy clustering, satellite fractions, cluster satellite luminosity functions, finding that 
subhalos reproduce well these observables. Our results imply that satellite subhalos 
corresponding to > 0.2 galaxies must be well-resolved down to 1 — 3% of their mass 
at infall to robustly trace the galaxy population. We also explore a simple analytic 
model based on dynamical friction for satellite galaxy infall, finding good agreement 
with our subhalo catalog and observations. 
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1 INTRODUCTION 

In the standard picture of galaxy formation, galaxies 
form at the centers of dark matter halos as baryons cool 
and contract toward the minimum of a halo's potential 
well dWhite fe Rees!ll978l : iBlumenthal et al.lll98l iDubinskil 
Il994l : |Mo et alJll998h T~But halos are not isolated objects: 
they merge over time, and smaller halos can survive as 
substructure halos (subhalos) of large r halos after infall 
dGhigna et al.lll998l; iTormen et all 1 19981 : iMoore et al"1ll999l ; 
iKlvpi n et al.l 1 19991 1. Thus, while galaxies form within dis- 
tinct host halos at high redshift, as the Universe evolves 
galaxies then correspond directly to dark matter subhalos 
(and since galaxies can be centrals or satellites, in our ter- 
minology a "subhalo" refers to a satellite substructure or 
the central halo). 

The evolution of subhalos involves more complicated 
dynamics than the formation of halos themselves. Satellite 
subhalos experience severe mass stripping as they orbit, and 
dyn amical friction causes their obits to sink toward halo cen- 
ter ()Ostriker fe Tremaindll975h . The evolution of galaxies is 
even more complex, involving gas dynamics, star formation, 
and various forms of feedback. Thus, the precise evolution- 
ary relation between galaxies and subhalos remains uncer- 



tain. In particular, what is the relation between galaxy stel- 
lar mass and subhalo dark mass? Do satellite galaxies expe- 
rience appreciable star formation after infall? How long after 
infall does a satellite galaxy survive, and what defines its fi- 
nal fate? How much of its stellar mass is funneled into the 
central galaxy in a merger, as opposed to tidally disrupting 
into Intra-Cluster Light (ICL)? 

Various prescriptions exist to map galaxies onto dark 
matter subhalos. The simplest is based on subhalo abun- 
dance matching (SHAM), which assumes a monotonic rela- 
tion between subhalo mass or circular velocity and galaxy 
stellar mass, populating subhalos such tha t one reproduces 
the observed stellar mass function (SMF: IVale fe Ostrikerl 
120061 : IConrov et al]|2006|; IShankar et al.ll2006l) . Alternately, 
semi-analytic models track the star formation histories of 
subhalos across t ime with em pirically motivated analytic 
prescriptions (seefBaugh 2006, for a recent review). These 
methods all involve assumptions regarding satellite galaxy 
star formation after infall, the correspondence of galaxies 
with subhalos in the case of severe mass stripping, and the 
eventual fates of satellite galaxies. 

Most models assume that all satellite galaxies eventu- 
ally merge with their central galaxy. Indeed, infalling satel- 
lites have long been thought to affect the mass and morpho- 
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logic al evolution of central galaxies l|Hausman fc Ostrikerl 
Il978h . However, galaxy clusters are observed to contain an 
appreciable amount of stellar material in the form of ICL, 
though given its low surface brightness, constraints on the 
ICL vary considerably, from 5 — 50% of the total cluster 
light (|Lin fe Mohrll2004l; [zibetti et al.l l 2005l ; iGonzalez et al.l 
120071 : Krick fc Bernstein! l2007l) . Additionally. ICL has been 
observed to contain signifi cant structure, including tidal 
streams l|Mihos et al.l l2005h . Thus, some satellite galaxies 
are at least partially, or perhaps entirely, disrupted into the 
ICL. 

In simulations, the criteria for subhalo merg- 
ing/disruption are influenced by mass and force resolu- 
tion. Insufficiently resolved subhalos will disrupt artifi- 
cially quickly , lead ing to the problem of "over-merging" 
ijKlvpin et al.l li999). A number of authors have examined 
in de tail mass stripping from satellite subhalos over time 
(e.g.. iMacieiewski et al.i r2009"l. however, using a fixed mass 
resolution limit for all subhalos, regardless of their infall 
mass, leads to lower mass subhalos becoming (artificially) 
disrupted more quickly. Cosmological A-body simulations 
are now attaining sufficient mass resolution to track sub- 
halos through many orbits and extreme mass stripping, in 
some cases to masses much smaller than the luminous mass 
of the galaxies they would host. At infall, the dark mass 
of a subhalo is an order of magnitude higher than its stel- 
lar mass, so tidal stripping of mass from subhalos cannot 
correspond directly to stellar mass stripping of its galaxy. 
Thus, in using dark matter subhalos to track galaxies, one 
must be careful both to treat subhalo merging/disruption in 
a self-consistent manner not dependent merely on numerical 
resolution, and to calibrate the point at which the galaxies 
they host become merged/disrupted. 

A variety of schemes have been used to define satellite 
galaxy merging/disruption via subhalos in simulations. Un- 
der the assu mption that a galaxie s survives as long as its 
subhalo does, IWetzel et all <|2009al lbh use an absolute mass 
threshold for all sub halos, corresponding t o their resolu- 
tion limit, similar to iKravtsov et ah! ((2004), whose model 
also implies that subhalos coalesce if their centers come 
within 50 ft -1 kpc. Several analyses assume that a satellite 
galaxy survives the disruption of its subhalo and eventu- 
ally merges with its central galaxy by imposing an ana- 
lytic infall timescale to a galaxy after subhalo disruption 
ifepringel et al.ll200ll: iKitzbichler fc Whitell2008l: ISaro et al.l 



120081 ; iMoster et al.ll2009l) . Alternately. IStewart et alj (|2009h 
allow a satellite galaxy to merge before its subhalo is dis- 
rupted. They assume that a galaxy's baryonic mass is cou- 
pled to the 10% most bound mass of its subhalo and define 
satellite merging when a satellite subhalo h as lost more than 
90% of its mass at infall. lYang et al . (2009) consider a model 
where the survival probability of a satellite subhalo is a (de- 
creasing) function of its infall mass to host halo mass ratio. 

Other criteria have been used in purely semi-analytica l 
mod els for the evolu t ion o f subhalos. IZentner et al.l (|2005l ) 
and iTavlor fc Babul l|2004l ) consider a satellite subhalo to 
be tidally disrupted when the mass of the subhalo has been 
stripped to a value less th an the mass within adi B r s , where 
r s is the subhalo's NFW dNavarro et al.l 1 1996h scale radius 
at infall, with aa la = 1 l|Zentner et al.ll2005l) or a^is = 0.1 
l|Tavlor fc Babul 2004]). Assuming a typical halo concentra- 
tion of c = 10, this implies that satellite subhalos will disrupt 



at 13% or 0.3% of their infall mass in the models, respec- 
tively - an order of ma gnitude difference. A s an alterna- 
tive to tidally disrupting, Zentner et alj (|2005l) consider the 
satellite subhalo merged with the cente r of its halo if it comes 
within 5 kpc. Applying the model of | Zentner et al.l (|2005f ) 
to the build-up of the ICL, l|Purcell et al.ll2007f ) allow satel- 
lite disruption to begin at ~ 20% of i nfall mass (see also 
iMonaco et~ai1l2006l : IConrov et aDl2007t for models of tidal 
disruption into the ICL). 

Finally, mass loss criteria have also been applied to 
modelling dwarf-spheroi dal galaxies in A-body simulations. 
iPefiarrubia et a l. (2008) find that dwarf satellites start to 
lose stellar mass to tidal stripping after 90% dark mass 
loss, but full disruption occurs only after > 99% dark mass 
loss. M accio et all (|2009h find good agreement between their 
semi-analytic model and the Milky Way satellite luminosity 
function using eidis = 0.5 — 1, corresponding to 5 — 13% of 
the infall mass remaining. 

Given the wide range in criteria for satellite merg- 
ing/disruption, it is informative to explore the implications 
of various criteria and to test empirically which match well 
to the observed galaxy population. Thus, in this paper we 
examine the extent to which galaxies can be mapped onto 
simulated dark matter subhalos throughout their evolution. 
We assign galaxy stellar masses and luminosities to subhalos 
by matching their number densities to observed galaxy sam- 
ples. We examine the effects that different criteria for merg- 
ing/disruption have on the satellite populations, including 
the HOD and radial distribution profile. By comparing spa- 
tial clustering, satellite fractions, and cluster satellite lumi- 
nosity functions to observations, we test the correspondence 
of galaxies to subhalos, and we constrain the criteria for 
satellite galaxy merging/disruption. 

Before continuing, we clarify the possible fates of satel- 
lite galaxies. The methods by which galaxies become re- 
moved from the satellite population are to (1) coalesce with 
the central galaxy, (2) coalesce with another satellite galaxy, 
(3) orbit outside their host halo, or (4) tidally disrupt into 
the ICL. As noted above, many galaxy evolution models as- 
sume only (1). As we describe below, our method for track- 
ing subhalos in simulations intrinsically incorporates (2) and 
(3). Thus, this paper focuses on methods for removing satel- 
lite galaxies via (1) and (4). Since our tracking of dark mat- 
ter subhalos does not unambiguously differentiate between 
these (and since both can occur simultaneously), we will use 
the term "removal" of satellite galaxies to indicate any mech- 
anism that causes (1) or (4), which we refer to as "merging" 
and "disruption" , respectively. 



2 NUMERICAL METHODS 

2.1 Simulations and Subhalo Tracking 

To find and track halos and their subhalos, we use dark 
matter - only A -body simulations using the TreePM code of 
IWhitel (|2002h . To track subhalos with high resolution, we 
use a simulation of 1500 3 particles in a periodic cube with 
side lengths 200 h' 1 Mpc. For our ACDM cosmology (Q m = 
0.25, fi A = 0.75, h = 0.72, n = 0.97 and ag = 0.8), in agree- 
ment with a w i de ar r ay of observations (ISmoot et al. 1992t 
iTegmark et al ] l2006t iReichardt et~aH [2009; Dunklev et al.l 
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2009), this results in particle masses of 1.64 x 10 8 /i _1 M 
and a Plummer equivalent smoothing of 3/i _1 kpc. Initial 
conditions were generated by displacing particles from a reg- 
ular grid using 2LPT at z = 250 where the RMS is 20% of 
the mean inter-particle spacing. We stored 45 outputs evenly 
in ln(a) from z = 10 to 0, with an output time spacing of 
~ 600 Myr at z — 0.1, the main redshift of interest in this 
paper. Given the limited volume of this simulation, and its 
implications for spatial clustering, we also use a separate 
simulation of 1500 3 particles of size 720 h^ 1 Mpc and a par- 
ticle mass of 7.67 x 10 9 H~ 1 Mq with the same cosmology 
and halo finder, though independent initial conditions. We 
use this larger simulation to populate our subhalo catalog 
onto a more accurate halo mass function with more accu- 
rate large-scale clustering. To examine the influence of cos- 
mological parameters, we also populate our subhalo catalog 
into a simulation of 1024 3 particles of size 500 ft" 1 Mpc with 
n = 1.0 and as ~ 0.9. We explore the influence of simulation 
size and cosmology in the Appendix. 

Our sub halo finding and track ing details are discussed 
extensively in lWetzel et alJ lf2009a.fi. though we have imple- 
mented a number of improvements for better performance, 
as described below. We find subhalos by first generating 
a catal og of halos using the Friends-of-Friends (FoF) algo- 
rithm l|Pavis et al.lll985T l with a linking length of b = 0.168 
times the mean inter-particle spacing|j We keep all groups 
that have more than 50 particles, and halo masses quoted 
below are these FoF masses. Within these "(host) halos" we 
then identify "subhalos" as gravitationally self-bound aggre- 
gations of at least 50 particles bounded by a density saddle 
point, using a new im plementation of the SubGnd algorithm 
l|Springel et alJlioOlh . where densities are smoothed over 32 
nearest neighbors. The "central" subhalo is defined as the 
most massive subhalo within its host halo, and it includes 
all halo matter not assigned to "satellite" subhalos, which 
is typically ~ 90% of the mass of the host halo. Thus, every 
sufficiently bound halo hosts one central subhalo and can 
host multiple satellite subhalos. A subhalo position is that 
of its most bound particle, and we define the halo center as 
that of its central subhalo. 

Each subhalo is given a unique child at a later time, 
based on its 20 most bound particles. We track subhalo 
histories across four consecutive outputs at a time since 
subhalos can briefly disappear during close passage with or 
through anot her subhalo. I n an im provement to the track- 
ing method of lWetzel et al.l (|2009aT ). we now interpolate po- 
sitions, velocities, and bound masses of these temporarily 
disappearing subhalos so that our tracking scheme does not 
artificially underestimate the satellite populat ion. 

Also, as discussed in IWetzel et all (|2009al ), subhalo his- 
tories often include cases of "switches" : if a satellite becomes 
more massive than the central, it becomes the central while 
the central becomes a satellite. In these cases, there is typi- 
cally not a well-defined single peak that represents the center 
of the halo profile, and the nominal center can switch back 



1 While of similar mass resolution, the sim ulation here has large r 
volume and better force resolution than in IWetzel et all ll2Q09ah . 

2 The longer linking length of b = 0.2 is often used, but it is more 
susceptible to joining together distinct, unbound structures and 
assigning a halo that transiently passes by another as a subhalo. 



and forth across output times. To avoid ambiguity in as- 
signing subhalo bound and infall mass, our tracking method 
now eliminates switches by requiring that once a subhalo is a 
central, it remains so until falling into a more massive halo. 
This history-based approach assigns as central the oldest 
subhalo in the host halo, which would typically correspond 
to the most massive and evolved galaxy, better mapping to 
how central galaxies are defined observationally. However, 
since this criterion is based on subhalo history, while Sub- 
find only uses the instantaneous density field to define the 
central (most massive) subhalo, the two definitions of sub- 
halo centrality are in conflict in a small fraction cases. In 
these cases, we defer to the history-based approach, and we 
swap the subhalo bound masses of the history-based central 
and Subfind-based central, such that the history-based cen- 
tral is always the most massive subhalo, and by extension 
the most massive galaxy. We find this produces more stable 
subhalo catalogs, in terms of individual mass histories and 
the resultant subhalo mass function. 



2.2 The Stellar Mass of Subhalos 

Since our goal is to test the correspondence between sub- 
halos and galaxies, a critical issue is how to relate stellar 
mass to a subhalo. The simplest prescription assigns stellar 
mass to a subhalo based on its instantaneous bound dark 
mass or maximum circular velocity. However, this method 
leads to satellite subhalos with significantly less concen- 
trated radial distribution s than seen for satell ite galaxies 
around the Milky Way dKravtsov et all l2004t) or galaxy 
clusters (|Piemand et al J 12004 ). Better agreement with ob- 
servations has been found using semi-anal ytic models of 
star f ormation applied to subhalo catalogs (|Springel et al.l 
l200ll ; IPiaferio et al.ll20o"lh . which correlate stellar mass to 
a subhalo's mass at the time the stellar component was as- 
sembled (before infall), thus indicating that it necessary to 
use subhalo mas s histories to accurately assign stellar mass 
[|Gao et al.ll2004h . 

The motivation for this approach comes from the 
following physical model. As a halo falls into a larger 
one and becomes a subhalo, tidal stripping quickly re- 
moves bound mass from the outer regions of the sub- 
halo, giving rise to significant mass loss shortly after ac- 
cretion. This behavior has been explored both in sim- 
ulations (e.g.. | Diemand et al. 12007 1_ Limousin et al.l 120091 ; 



iPenarrubia et al.l 120081; IWetzel et al.l l2009ah and obser- 
vatio n ally through ga l axy-g a laxy lensing ([Limousin et al.l 
120071 ; iNataraian et al.l 120071 l2009h . However, the galaxy, 
being more compact and at the center of the subhalo, re- 
mains intact longer. Thus, dark mass stripping does not di- 
rectly correspond to stellar mass stripping. The subhalo's 
gaseous component - being more diffuse and readily affected 
by ram-pressure stripping, tidal shocks, and adiabatic heat- 
ing in the hot group/cluster environment - quickly becomes 
stripped after infall, halting subsequent star formation in 
the satellite galaxy. Thus, galaxy stellar mass is expected to 
correlate with the subha lo's mass immediately before infall 
|Nagai fc Kravtsovll2005h . 

Observationally, there are a variety of claims for the 
star formation efficiency of satellite galaxies, with signif- 
icant systematic dependence on group/cluster definition 
or HOP parametrization, criteria for blue/red (star form- 
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ing/ quenched) split, and inclination effects fsee lSkibball2009l . 
for a comparison of HOD color models). The HOD analy- 
sis of the Sloan Digital Sky Survey (SDSS) of IZehavi et al.l 
(2005) implies few faint blue satellites, and that at higher 
luminosity blue satellites are 3 — 5 times less common than 
red sa tellites in groups . Exa mining group catalogs from 
SDSS, IWeinmann et all (|2009h find an appreciable fraction 
(> 20%) of satellites in groups are blue, and fit a simple 
fading model where star fo rmation in satellites quenche s 
~ 2 Gyr after infall (see also lKang fc varTde n Bosch 2008). 
Using a subhalo catalog wit h an analyti c pres cription for 
star formation and quenching, IWang et all (|2007h find an ex- 
ponential decay time for satellite star formation after infall 
of ~ 2 .3 Gyr. At z ~ 1, the clustering results of ICoil et al.l 
|2008l ) imply that immediate star formation quenching at 
infall is likely too restrictive at higher redshift. 

Hydrodynamic simulations also indicate short yet finite 
timescales for halting star formation. High-resolution SPH 
simulations at z ~ indicate that sate llite star formatio n 
can continue for up to 1 Gyr after infall (ISimha et alj [2009). 
While the assumption that satellite gas accretion rates are 
significantly lowe r than those of ce ntral galaxies may break 
down at z > 1 l|Keres et al 1 l2009h . the significantly more 
recent infall times of s atellites at high red s hift m ediates this 
concern. Promisingly, iNagai fc Kravtsovl (|2005h find good 
agreement in comparing the radial profiles of subhalos in 
dark matter-only simulations selected on infall mass with 
galaxies selected on stellar mass in dark matter plus hydro- 
dynamic simulations. 

Despite the above uncertainties, attempts to use sub- 
halo infall mass (or infall maximum circular velocity) to 
assign stellar mass have been successful at reproducing 
various observed ga l axy properties at a r ange of redshifts 
llBerrier et"aH 120061: IConrov et all l2006t IVale fc Ostrikerl 
l200d : IWang et al.ll2006l ; lYang et al.ll2009h and so this is the 
approach we use here. The success of this method may be a 
derivative of its implicit allowance for some star formation 
after infall, as we describe below. Furthermore, our analy- 
sis is restricted to z < 1, where the approximation of rapid 
quenching of star formation is most valid. 

2.3 Assigning Infall Mass and Stellar Mass 

Motivated by its expected correlation with stellar mass, our 
scheme for assigning subhalo infall mass, M- m f, is as fol- 
lows. Centrals subhalos are assigned Minf as their current 
bound mass. This mass represents halo mass not bound 
to any satellite subhalos and so it should track the gas 
mass available to accrete onto the central galaxy. We as- 
sign to a satellite subhalo its mass when it fell into its cur- 
rent host halo (its subhalo mass when it was last a cen- 
tral subhalo). With the expectation that a subhalo merger 
corresponds to a galaxy merger, if two satellite subhalos 
merge, their satellite child is given the sum of their infall 
masses. Thus, our tracking scheme inherently incorporates 
stellar mass growth though satellite-satellite mergers, an 
often -ignored but non-tr ivi al contribution to ga laxy evolu- 
tion I Wetzel et al.ll2009allbl;[Seek Kim etal]|2009l ). If a satel- 
lite subhalo's orbit brings it outside its host halo, thus be- 
coming its own distinct smaller halo, the subhalo is assigned 
its satellite parent's infall mass during the output immedi- 
ately following ejection, but is assigned its current bound 
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Figure 1. Relation between stellar mass, M star , and subhalo in- 
fall mass, Msui^inf, using subhalo abu ndance matching against 
the observed stellar mass function of ICole et all (|200ll ). Solid 
curve shows relation to our full subhalo catalog (no removal 
threshold), dotted curve shows the relation under our maximal 
threshold for satellite removal, and dashed curve shows the mean 
relation when using a log-normal 0.2 dex scatter in M s t ar at a 
fixed Msub^inf. Inset shows ratio M s tar/M su b,inf . 



mass for subsequent outputs if it remains a separate central 
subhalo. 

We assign stellar mass to our subhalo catalog using sub- 
halo abundance matching (SHAM) such that by construc- 
tion we recover the observed stellar mass function. We rank 
order our subhalo catalog by infall mass, and assign stellar 
mass such that n(> Mj n f) = n(> M a t ar ), u sing the observed 
stellar mass function of ICole et all (j200ll ). Fi gure [1] shows 
the M s tar — Mi n f relation, for our full satellite subhalo cat- 
alog (no removal threshold) and for the maximal threshold 
for satellite removal we consider. Since the satellite fraction 
increases at lower mass, the discrepancy is strongest there, 
though the offset is co mparable to the error in the observed 
stellar mass function (jCole et al.ll200ltl . The star formation 
efficiencies in Fig. [T] are in good agreement with the con- 
strains from weak lensing of iMandelbaum et al.l (2006) and 
also agree with the Milky Way, which has stellar mass of 
5 x W 10 Mq and a dark halo mass of about 2 x 10 12 M n 



llBinnev fc Tremainel [20081) . making M sta r/M ha i « 3%. See 
iMoster et al.l 1 20091 ) for a detailed comparison of derived 
Mgtar — M au bhaio relations in the literature. 

Despite the seeming requirement for rapid satellite 
quenching after infall, SHAM does implicitly allow for satel- 
lite star formation. This is because we match stellar mass 
to satellites' infall mass at an observational epoch, not at 
the time of infall. At z < 1, satellite times since infall are 
typically several Gyr, and during that time satellite M- ln f re- 
mains fixed. Thus, if there is appreciable evolution (growth) 
in the global M star — M- m { relation during that time, this 
will manifest itself as stellar mass growth for satellites. We 
have explored additional satellite star formation after infall 
by increasing all satellite infall masses by a fixed fraction, 
but we find that this is simply degenerate with increasing 
our satellite removal threshold given in the next section. 

While the SHAM technique assumes a one-to-one cor- 
respondence between Mj nl and M star , there may be consid- 
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erable scatter in the relation. Fits to the spatial clustering 
of luminosity or stellar mass limited galaxies samples fre- 
quently require a "soft" turn on in the central galaxy oc- 
cupation least they overpredict the large-scale clustering, 
which would naturally arise from scatter in observab le at 
fixed (sub)halo mass. For example. IZheng et al.l (|2007l ). ap- 
plying Halo Occupation Distribution (HOD) modeling to 
galaxies at z ~ and ~ 1, find log-normal scatter in lumi- 
nosity at fixed halo mass of 0.15 — 0.3 dex, decreasing with 
increasing halo mass, consistent with the scatter of ~ .15 
dex found in the group catalogs of | Yang et al.1 (|2008h at 
z ~ 0. Similarly, Ivan den Bosch et al.l l|2007f ) find a lower 
limit of 0.2 dex scatter in halo mass at fixed luminosity from 
conditional lu minosity fu n ction modeling, consistent with 
the results of [ More et all (|2009f) from satellite kinematics. 
iTasitsiomi et al.l l|2004T l find good agreement in spatial clus- 
tering when matching subhalos to galaxies if they imposed 
significant scatter (0.6 dex) in luminosity at fixed subhalo 
maximum circular velocity. 

To examine the importance of this scatter when com- 
paring against observed galaxy samples, we also produce a 
stellar mass catalog by imposing a fixed 0.2 dex (log-normal) 
scatter in stellar mass at fixed subhalo infall mass such that 
we recover the observed SMF. The mean M sta r at a fixed 
Minf is shown as the dashed curve of Fig. [T] Since the mass 
function steeply falls with mass, introducing scatter biases 
M s t a r to a lower value at a fixed Minf. Higher scatter also 
causes a stellar mass threshold to correspond to a lower ef- 
fective subhalo Mi n f threshold This effect is particularly 
strong above Mint « 10 12 /i _1 Mq, where the subhalo mass 
function starts to fall exponentially with mass. 

We follow a similar approach to SHAM when we com- 
pare to observations of magnitude limited samples (Q, 
matching the number densities as n(> Minf) = n (> L). 
The two assumptions - that stellar mass and luminosity are 
set by infall mass - are not fully consistent since luminos- 
ity may evolve even if stellar mass does not. However, it is 
not unnatural to assume the rank ordering is preserved in a 
population sense. 



3 IMPACT OF SATELLITE REMOVAL 

As a subhalo orbits in its host halo, dynamical friction re- 
moves energy from its orbit, bringing it closer to the central 
galaxy, while tidal forces strip its mass. While the satellite 
galaxy is expected to retain its stellar mass throughout most 
of the subhalo's mass stripping, eventually (though possi- 
bly in longer than a Hubble time) the galaxy will become 
removed from the satellite population, either from merg- 
ing with the central galaxy or becoming disrupted into the 
diffuse ICL (or both). We now examine the impact differ- 
ent removal criteria have on the satellite population. We 
first focus on physicall y informative properties in the con- 
text of the halo model liseliakll2000l: IPeacock fc Smithll200d : 
iBerlind fc Weinberdl2002l ; ICoorav fc Shethll2002l) . including 



3 This is not immediately apparent in Fig. where it might 
appear that introducing scatter causes increased effective Mi n f 
at fixed Mstar- This is because we show the mean M s t ar at a 
given Mi n f, which is qualitatively different than the mean Mi n f 
at a given Mstar ■ 



the HOD and radial profile. However, since these are not 
direct observables, we reserve detailed comparisons with ob- 
servations, via spatial clustering and satellite fractions, to 
the next section. 

As described in the introduction, various prescriptions 
for satellite subhalo removal have been used. We focus pri- 
marily on the one used most often in previous work, which 
is also the most straightforward given our tracking: remove 
subhalos where the instantaneous bound mass to infall mass 
ratio, /inf = Mbound/Minf , falls below some threshold. While 
dark matter-only simulations do not incorporate the com- 
plex hydrodynamics involved in galaxy evolution, if a satel- 
lite galaxy's gas is stripped before its stellar component, 
then it is conceivable that a simple criterion on subhalo dark 
matter stripping yields a good approximation for galaxy 
stellar mass removal since both experience purely gravita- 
tional interactions. We examine a range in threshold /i n f 
from 0.01 — 0.1, which as we show in ^brackets the obser- 
vations, and even subhalos stripped to /i n f = 0.01 remain 
well-resolved in the mass ranges we consider below. We re- 
fer to our subhalo catalog with no removal threshold as the 
"full" tracking model. Based on convergence tests, our sub- 
halo catalog at z — 0.1 is robust to artificial disruption down 
to ~ 10 11 ' 5 H~ 1 Mq, so the full tracking model serves as an 
upper limit to the satellite population. We explore different 
removal criteria in SjfJ] 

3.1 Halo Occupation Distribution 

Figure [5] (top) shows the Halo Occupation Distribution 
(HOD) for our full tracking model ("no /inf") and for dif- 
ferent minimum /inf = Mb OU nd/Mi n f thresholds. We fit the 
HOD to 



N m 



( M ha 



V Mi 



-M cut /M halo 



(1) 



with Mi, M cut and 7 as free parameters. The slope, 7, is 
given in Fig. [2] for each value of /i n f. Raising the threshold 
for satellite removal, /i n f, affects satellites more in higher 
mass halos, leading to a shallower slope. As the threshold 
increases, the HOD shoulder also becomes broader, mean- 
ing that more halos host only one galaxy. This indicates that 
satellite subhalos of a given Minf are stripped more in higher 
mass halos, because of the stronger tidal fields and longer 
dwell times. This implies that the details of modeling satel- 
lite galaxy removal are more critical in cluster-mass halos 
than in lower mass groups. 

The dotted curve also shows the HOD if satellites are 
never removed by disruption or merging with the central 
galaxy, in other words, every infalling galaxy survives until 
z = 0.1 (which we measure by analytically retaining ev- 
ery infalling halo). This provides a strict upper limit to the 
HOD, contingent only upon the accuracy of the halo merger 
trees and not on subhalo tracking and resolution. As we ex- 
plore in S|4] this scenario is highly disfavored. 

In the halo model, the spatial clustering on large scales 
is dominated by galaxy pairs in separate halos ("2-halo" 
term), which depends on the first moment of the HOD, 
(A)(Mhaio), where N = Neat + 1. However, the cluster- 
ing on small scales is dominated by pairs within a halo 
("1-halo" term), which depends on the second moment, 
(AT(JV-l))(M ha io). Figure E] (bottom) shows the normalized 
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Figure 2. Top: Halo Occupation Distribution (HOD) for satel- 
lite subhalos with Af inf > W 12 h- 1 M Q (M star > 10 10 ' 8 M Q , 
L > 0.6 L t ), for different minimum removal thresholds, /; nf = 
A^bound/A^inf ■ Solid curve shows all resolved satellites, regard- 
less of mass stripping ("no /; n f"), while dotted curve shows all 
satellites ever accreted, even if they fall below the resolution limit 
(a model in which satellites never merge or disrupt). 7 indicates 
best-fit slope to Eq. ^ Errors are standard deviation of the mean 
HOD at fixed halo mass. Bottom: Normalized second moment 
of satellite HOD. 



second moment of the satellite HOD, denned as (JV 8a t (JV Ba t — 
V)} / (N S!i t) 2 . The normalized second moment asymptotes to 
~ 1 at high halo mass, reflecting a nearly Poisson distri- 
bution, and it approaches zero at low halo mass where the 
probability of hosting one satellite remains finite while the 
probability of hosting two satellites goes to zero. Figure [5] 
shows some dependence on removal threshold, particularly 
at high halo mass where raising the threshold increases the 
second moment. This is driven by correlated infall, such as 
multiple satellites infalling as a group, and thus correlated 
amounts of mass stripping. This is supported by the fact 
that the "no removal" scenario (dotted curve), which rep- 
resents all satellites ever accreted and is insensitive to mass 
stripping, is closer to Poisson. 
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Figure 3. Top: Satellite radial density profile for satellites with 
M^t > 10 12 h- x M @ (M s tar > 1O 1O - 8 M , L > 0.6L,). Solid 
curve shows all resolved satellites (no removal threshold), while 
subsequent curves show different minimum removal thresholds. 
Radii are defined by distance to central galaxy and are scaled 
to the host halo virial radius, R v i r . Densities are normalized to 
the number within 2 i? v i r using no threshold. Errors bars indicate 
Poisson error in each radial bin. Bottom: Ratio of profile using 
a given /; n f threshold to that of using no removal threshold. 



3.2 Satellite Radial Density Profile 

In addition to the number of satellites in a halo, we also 
explore how removal thresholds affect their locations. Fig- 
ure|3](top) shows the satellite radial density profile, obtained 
from stacking all halos hosting the given satellite subhalos. 
Radii are defined by the distance of a satellite to its central 
galaxy and are scaled to the host halo's virial radius, R v i T , 
obtained from the halo FoF mass and concentration assum- 
ing a spherical NFW density profile (R v i r ~ R200C at z m 0). 
Since halos are triaxial, particularly for those with signifi- 
cant substructure (recent mergers) , the assump tion of spher- 
ical symmetry is highly approximate (see e.g.. IWhitdl2002l . 
Fig. 2). Additionally, subhalos preferentially lie in a plane 
defin e d by their host halo's major axis (e.g., Kroupa et al.l 



— a.T— 

120051 : iFaltenbacher et al.l 120081 : iLibeskind et all 20091) . be 
havior also found for galaxies in SPSS terainerdll2005h and 
dwarf satellites around the Milky Way l|Metz et alj|20615h . 
This implies that the assumption of spherical symmetry is 
even less valid for subhalo populations, and that a spheri- 
cal overdensity halo finder will miss satellites aligned along 
the halo major axis. Because of these effects, we find that 
selecting satellites out to _R v i r only captures 70% of the en- 
tire satellite population, while going out to 2 R v i T captures 
95%. Hence, number densities are normalized by the counts 
(using no removal threshold) within 2 i? v ir • 

The full subhalo tracking model (no removal threshold) 
shows satellite subhalos tracing the NFW profile to good 
approximation down to 0.01i? v ir, though this may repre- 
sent over-resolution. Since stripping is more prominent in 
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the halo's dense central region, raising the threshold prefer- 
entially removes satellites from the halo center, leading to 
a less concentrated profile. Figure [3] (bottom) shows the ra- 
tios of the profiles using various /j n f values to that of the 
full tracking model. In all cases, significant suppression oc- 
curs below 0.4-R v i r . Interestingly, for no /i n f threshold does 
the profile asymptote to unity, implying that mass stripping 
affects satellites even at large radii. The asymptotic value 
has no significant dependence on Mj n f or on time since in- 
fall, indicating that this effect is not driven by satellites that 
have passed through the halo center, but instead by subhalos 
being significantly stripped during or soon after infall. 

These profiles can be compared with observation- 
ally measured satellite surface density profiles in galaxy 
groups/clusters, a method typically based on stacking the 
profiles of galaxies from several groups/clusters. Some work 
indicates that satellite galaxies trace well their host halo's 
NFW profile. Using a compilation of local galaxy clusters, 
iDiemand et al.l l|2004f l find that the surface density profiles 
of their member galaxies match that cluster-mass halo parti- 
cles in simulations within 20% s catter and no systematic off - 
set down to 0.05 R v i T . Similarly. Ivan den Bosch et al.l (|2005h 
examine satellite galaxies in Two Degree Field Galaxy Red- 
shift Survey (2dFGRS), finding that they are consistent with 
following their host halo's density profile, though they were 
not able to discriminate distributions in halo cores. 

Other work indicates that, while satellite density pro- 
files can be fit by an NFW profile, they require a lower con- 
centration parameter than the dark matte r. Using clusters 
of mass > 3x I0 13 fo^M^- fLm et alj (|2004l) find that galaxy 
profiles down to 0.02 R V1I require an NFW con centration pa- 
ramet er ofc = 2.9±0.2. A similar analysis bv lMuzzin et al.l 
120071 ) finds that galaxies trace NFW with c = 4.1 ± 0.6 
down to 0.5-Rvir- Examining lo wer mass (< l O ^ft-^ Mfl) 
galaxy groups/clusters in SDSS, lHansen et al.1 (|2005h find 
that satellite profiles down to 0.1-R v ir can be well- fit with 
NFW profiles using c < 3, with sim ilar results from the 
2dFGRS and SDSS group catalogs of l|Yang et al.ll2005l l. 

Concentrations of c = 3 — 4 are lower than the average 
halo concentration of c = 5 — 6 for cluster-mass halos in 
our simulations. Interestingly, an NFW profile with c = 3 is 
consistent with the satellite distribution in our simulation 
down to 0.06 i?vir using /i n f = 0.01. However, differences in 
halo virial radii estimates, density normalizations, cluster 
centering, and galaxy luminosity thresholds make a detailed 
comparison difficult, as highlighted by the diversity in con- 
centrations derived observationally. 

One strong systematic uncertainty in measuring the 
profiles is how one chooses the halo center. Two common 
methods of defining a cluster center are the location of the 
Brightest Cluster Galaxy (BCG) and the peak of X-ray 
temperature profile. However, these are often offset by sev- 
eral 10's of kpc. For example, for rich clusters in the SDSS 
MaxBCG catalog, the me dian offset between the BCG and 
X-ray center is 58 ft -1 kpc l|Koester et al.|[2007f ). To estimate 
how such uncertainties propagate, we examine the degree to 
which centering changes the measured density at 0.1 7i v ir by 
giving the central subhalo an arbitrary offset. We find that 
an offset of 10, 25, or 50 h~ l kpc leads to a change of +1%, 
+4% and —4% in the measured density, the error growing 
rapidly with offset amplitude. Requiring less than a 10% re- 
duction in measured density at 0.1 i? v ir for clusters requires 
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Figure 4. Satellite radial distribution probability (not divided by 
volume) at the time of crossing below a given removal threshold, 
/ inf , for satellites with Af inf > 10 12 fe _1 M Q (Mstar > 10 10 ' 8 M Q , 
L > 0.6 L*). For reference, solid curve shows radial distribution 
of all resolved satellites at z = 0.1. 

an offset of less than 60 h^ 1 kpc in the central subhalo. Such 
offset may be minimized by selec ting BCGs based on their 
color profiles (jBildfell et al.ll2008h . 

3.3 Radius at Removal: Merger vs. Disruption 

While Fig. [3] shows the profile of extant satellites, it does 
not indicate the location of satellites when they became 
removed and where their central particles, associated with 
a stellar population, end up. This would provide a mea- 
sure of whether satellite galaxies merge with the central 
galaxy or disrupt into the ICL. One route to examine this 
would be to track central "star" particles in subhalos after 
they become re moved and examine t heir overall distribu- 
tion at z ~ dMurante et all 120041: IWillman et al.1 1 20041 : 
ISommer-Larsen et al.ll2005l ; iRudick et al.ll2006f) . We take a 
different approach and examine galaxy radial distribution at 
the time of removal. While Fig. [3] shows satellites above a 
threshold /j n f, we now examine satellites as they fall below 
the threshold by selecting those within 25% of a given /i n f 
value. 

Figure U shows the probability distribution that a satel- 
lite is at a given scaled radius for a given /i n f value. For 
reference, the solid curve shows the radial distribution of all 
resolved satellites at z — 0.1, regardless of stripping. These 
are obtained from stacking all halos hosting the given satel- 
lites, though the results do not change significantly if we look 
only at cluster-mass halos. Satellites fall below a lower /i n f 
threshold at a smaller radius. This is expected in a model 
where subhaloes are stripped of mass as dynamical friction 
brings their orbits to smaller radii. All thresholds display 
a broad distribution with radius, since stripping will also 
depend on satellite orbital p arameters and time since infall 
l|Bovlan-Kolchin fc Mai 120071 ). However, the distribution is 
more peaked for higher /i n f thresholds. This likely arises be- 
cause significantly stripped subhalos are closer to halo cen- 
ters, where their velocities are both higher and more radial. 
Indeed, as the removal threshold is raised from /i n f = 0.01 
to 0.1, the average (absolute) ratio of radial to tangential ve- 
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locities for the satellites drops from 1.5 to 1.1. The average 
for the entire satellite population is between these at 1.4. 

To understand their fates, one must know the direc- 
tion in which satellites are moving as they are removed. If a 
satellite disrupts at large radius but is moving toward halo 
center on a highly radial orbit, its stellar mass may still fun- 
nel to the central galaxy. As a measure of this, we examine 
the fraction of satellites that are moving net inward towards 
halo center as they cross below /i n f. We find only mild de- 
pendence on /i n f: as the removal threshold is raised from 
/inf = 0.01 to 0.1, the fraction of satellites moving inward 
falls from 58% to 52%. Both are lower than the average for 
the entire satellite population of 62%. 

The above results imply that more highly stripped satel- 
lites are more likely moving inward and their orbits are more 
radial. However, both of these trends with mass stripping are 
mild. Fig. [4] shows that, even in our most conservative case 
of /inf = 0.01, the (broad) peak of the radial probability 
distribution occurs at ~ 50 h~ kpc. This, coupled with the 
fact that almost half of satellites at removal are moving out- 
ward and have significant tangential velocity components, 
implies that a significant fraction of satellite galaxies dis- 
rupt into the ICL, and d o not immediately merge with the 
central galaxy (see also iMonaco et al.| [20061 ; IConrov et al.l 
120071 ; iPurcell et al.ll2007t IWhite et alj|2007h . 



100 r- 



3.4 Analytical Model for Satellite Removal 

As an alternative to numerically tracking satellite subhalos, 
we can assign infalling halos a removal time at accretion 
and examine what impact different timescales have on the 
HOD. While individual merging /disruption times depend 
on satellite orbital parameters dBinnev fc Tremaind 120081 : 
iBovlan-Kolchin et al.ll2008l : Ijiang et al.ll2008h . we use a sim- 
pler parametrization of satellite removal time which depends 
only on the halos' mass ratio at infall, assumed to hold for 
an ensemble average of satellites. Here, an infalling satellite 
halo with mass M sati i n f merges with the central galaxy or 
becomes disrupted on a timescale 



tdyn — Cdyn 



M h alo/M Ba t,inf 
ln(l + M ha l /M sati in 



^Hubble 



(2) 



where fHubbic = H~ 1 (z), Mhalo is the larger halo mass, and 
we leave Cd yn as a free parameter to match to different satel- 
lite removal thresholds. 

Figure [5] shows the HOD for multiple Mi n f-limited sam- 
ples and two removal thresholds, /j n f. Also plotted are the 
HODs using different values of Cd yn for dynamical friction. 
In each case, the Cd yn value is chosen to best match the sub- 
halo catalog using the given /i n f value. Impressively, this sole 
free parameter can match well the subhalo catalogs across 
a broad range of subhalo and halo masses, and for multi- 
ple removal thresholds. Additionally, since satellite removal 
times depend on their orbital parameters the relative ro- 
bustness of Cdyn across a wide halo and subhalo mass range 
indicates that satellite orbits do not vary significantly with 
mass, though Fig. [5] shows some hint of shorter infall times 
(more radial orbits) for satellites in more massive halos. As 
we shall see in Sj5] using Eq. [2] with the same value of Cdyn 
for the given /i n f also works well at higher redshift. 
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Figure 5. Halo Occupation Distributions (HOD) for 
multiple Afi n f -limited samples, using removal thresh- 
olds /inf > 0.01 (top) and / inf > 0.03 (bottom). 



Mint = 1O 11 ' 5 ,1O 12 - u ,1O 12 - 5 ,1O 13o /i- 1 M corresponds to 
M star = 10 10 - 3 ,10 10 - 8 ,10 ni ,10 113 (L = 0.2,0.6, 1.3, 2.1 L*). 
Also shown are HODs for different values of Cd yn f° r dynamical 
friction infall time, chosen to match the subhalo catalog. 



4 COMPARISONS WITH OBSERVATIONS 

Having examined trends in subhalo properties with varying 
removal threshold, we now seek to constrain this freedom 
by comparing with observations. While the HOD and ra- 
dial profile of the above sections are physically informative, 
they are not direct observables. The most direct comparison 
with data is via spatial clustering, which provides a scale- 
dependent test of our subhalo catalog at various masses. 
In addition, we compare to observed satellite fractions and 
cluster satellite luminosity functions, though these measures 
are less robust since they are derived from HOD modeling of 
spatial clustering or from constructing group catalogs, both 
of which are subject to systematic uncertainties in halo mass 
definition, HOD parametrization, and assumed cosmology. 

We again use subhalo abundance matching, as outlined 
i]2.3l but while we matched to the 2dF stellar mass function 
in i)2.3l in this section we match to the number densities 



of r-band luminosity threshold samples in SDSS in order 
to provide a robust comparison against observed clustering 
results. Thus, here the requirements on matching are less re- 
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strictive, in that we need only match to a single threshold at 
a time rather than reproduce the entire luminosity function. 

4.1 Spatial Clustering 

For spatial clustering measurements, we compare a gainst 
the SDSS galaxy clustering results of Ze havi et al.l (|2005f) 
at z » 0.1. To measure spatial clustering, we use the two- 
dimensional projected galaxy auto-correlation function 



w p (r p ) 



dn £(r p ,7r) 



(3) 



which has the advantage of effectively integrating over 
redshift-space distortions, which we incorporate using sub- 
halo velocities. We use 7r max = 40 /i -1 Mpc to match the 
measurements against which we compare. This value repre- 
sents a good balance between integrating over redshift-space 
distorti ons and minimizing correla ted noise from large-scale 
power |Padmanabhan et al.ll2007l) . 

Both large- and small-scale clustering are sensitive to 
the satellite population. In the limit that all galaxies are 
centrals (a halo hosts only one galaxy), the galaxy correla- 
tion function is simply that of the mass-limited halo sample 
of the same number density. As the threshold for satellite re- 
moval decreases and the satellite population increases, high- 
mass halos will host multiple galaxies. This serves to raise 
the minimum halo mass of the sample (at fixed number den- 
sity) while further weighting the clustering from high-mass 
halos which host more satellites. Satellites in a halo are also 
close spatially, to each other and to the central. Thus, satel- 
lites boost both the large- and small-scale clustering. 

The small-scale clustering is also dependent on the ra- 
dial profile of the satellites, though this sensitivity is mild. 
In the limit that central-satellite pairs dominate, the clus- 
tering on scale r p comes from satellites at a large fraction of 
r v i r in small halos and satellites nearer the center of larger 
halos. Given the steepness of the mass function and the avail- 
able halo volume, the first contribution dominates, weaken- 
ing the sensitivity to the inner profile. Therefore, w p (r p ) is 
only mildly sensitive to the radial dependence of satellite re- 
moval, though in the extreme limit of no satellites it drops 
off rapidly on small scales because of halo exclusion. 

Figure [6] shows w p (r p ) of luminosity-threshold samples 
from SDSS. These are compared against the Mi n f-threshold 
subhalo catalog for several removal thresholds, f- ln f, each 
matched to the threshold number density of the SDSS sam- 
ple. The magnitude-limited samples of M r < —20.5 and 
< —21.5 correspond to subhalo Mjnf -limited thresholds of 
Minf ~ 10 12 and 10 13 /i _1 M0, spanning a range we are 
able to probe with robust resolution and good statistics. 
Impressively, the simple prescription for satellite removal 
of /inf = 0.01 — 0.03, combined with simple number den- 
sity matching to assign luminosity, matches well at small 
and large scales, across over an order of magnitude in Mi n f . 
We find similar agreement with the M r < —20.0 sam- 
ple, correspo nding to M n f > 10 11 75 /i~ 1 Mp ) , if the "Sloan 
Great Wall" l|Gott et al J 120051 ) is removed, but because of 
the involved uncertainty we do not include this sample to 
constrain /i n f. The M r < —19.5 sample, corresponding to 
Mini > 10 11 ' 5 h~ 1 MQ, agrees as well but is on the edge of 
robust resolution, so it does not provide a robust constraint. 

We also investigate the influence of luminosity-mass 
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Figure 6. Projected auto-correlation function at z = 0.1 for 
several removal thresholds, /i n f, (curves) each matched in num- 
ber density to magnitude-limited samples in SDSS assuming 
no L r — Mj n f scat ter and compared with measured w p (r p ) of 
IZehav i et al. (20Q5|) (points). Top: w p (r p ) (top) and ratio of sim- 
ulation to observed w p (r p ) (bottom) for M r < —20.5 sample. 
Bottom: Ratio of simulation to observed w p (r p ) for M r < —21.0 
and < —21.5 samples. M r < —20.5, —21.0, and —21.5 correspond 
to subhalo M inf > 10 12 , 10 12 - 5 , and lO 13 h- 1 M . Also shown is 
the reduced \ 2 of the fit to observation for each /; n f . 



scatter. Because of the steepness in the subhalo mass func- 
tion, adding scatter lowers the effective Mi n f threshold at a 
fixed number density, and this effect is stronger at higher 
Minf (see Fig.HJ. Thus, Fig. [7]shows w p (r p ) for most lumi- 
nous samples, assuming 0.2 dex scatter in luminosity at fixed 
Mi n f. As expected from the decreased effective mass thresh- 
old, adding scatter reduces the clustering amplitude, and 
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Figure 7. Same as Fig. [6] (bottom), but luminosity is matched 
to subhalos assuming 0.2 dex scatter in L r at fixed M- ln {. 



the effect is stronger at farger scafes. However, the goodness 
of agreement between observations and f- m f — 0.01 — 0.03 is 
not significantly changed, except for the M r < —21.5 sample 
where agreement is improved somewhat, suggesting neces- 
sary scatter at high mass. 



4.2 Satellite Fraction 

We next compare against the observed satellite fraction: the 
fraction of galaxies above a given luminosity threshold that 
are satellites. Since the satellite fraction is an integral over 
the HOD, it integrates the effects of ij3]into a single number 
as a function of Mi n f. as the threshold for removal is raised, 
the satellite fraction drops. 

We compare against a variety of values derived from 
SDSS and 2dFGRS, all at a median z = 0.05 - 0.1. Again 
we stress that determining satellit e fractions from obs erva- 
tions requires sign ificant modeling. IZehavi et all (|2005h and 
IZheng et all {2007) satellite fractions were both derived from 
fitting HODs to the number density and spatial clustering 
of galaxy samples from SPSS, thoug h they used different 
HOD parametrizations. I Tinker et al.l l|2007l ) u sed a simi- 
lar tec hnique applied to 2dFGRS. Additionally, lYang et all 
(2008) satellite fractio ns were obt a ined directly from the 
SDSS group catalogs of lYang et ail (|2007f ). 

Figure [8] (top) shows observed satellite fractions, abun- 
dance matched to our subhalo catalog and plotted as a 
function of Mjnf. The observationally-based analyses yield 
similar results, which also agree with va lues from galaxy- 
galaxy lensing ijMandelbaum et al.l l2006h . though there is 
consi derable scatt e r. At the low mass end, the amplitude 
from lTinker et all (|2007T) is anomalously low, which may be 
driven by the fact that the 2dFGRS is a blue-selected sam- 
ples and thus could be missing a large population of satellites 
there, since satellites are pr eferentially redder than centrals 
l)van den Bosch et al.l l2008h . Given the consistency of the 



0.4 

c 

\ 

a 

£ 0.3 

c 
o 

4-1 

o 

£ 0.2 
t. 

m 



0.1 



(0 
CO 



-1 — I — I I I 1 1 1 



-1 — I — I I I 1 1 1 



-1 — I — I I I 1 1 1 



z=0.1 



no f,_ 



10" 



^ 0.4 

3 
C 



0.3 



c 
o 

'£ 
o 

10 

u 
fc. 



0.2 



0.1 



co- 
co 




_i i 



j i ' 



_i i 



10" 



10 la 



M, 



10 13 
[h-'M ] 



10 14 



Figure 8. Satellite fraction vs. infall mass. Top: Using various 
removal thresholds, /i n f, from the subhalo catalog. Points show 
observationally derived values, matched in numbe r density to 
our su bhalo c atalog, from | Ze havi et all j2005h fZ05). |Zheng etjil] 
(l2007h (ZOTL lTinker et al.l ||2007|) (T071, and lYang et al.1 ll2008t) 
(Y08). Bottom: Using the analytic model for satellite halo infall 
of Eq. [2] and various values of Cd yn for dynamical friction. Good 
agreement with observation (no rollover) continues at low mass 
since the analytic model depends only on resolving halos prior to 
infall. Dotted curve shows keeping all satellites ever accreted. 



three other samples at low mass, we consider their values 
more robust in th at regime. At the hi gh mass end, the satel- 
lite fractions from lZheng et all (|2007l ) show an upturn; since 
the satellite fraction is expected to decrease with mass and 
the amplitude of the highest two mass bins are considerably 
different than the other samples, we will not consider them 
as robust data points. 

Also shown are satellite fractions from our subhalo cat- 
alog for various removal thresholds, /i n f. Since the satellite 
fraction should rise with decreasing galaxy mass, the rollover 
at low mass clearly shows the limits of numerical resolution 
below M inf = 10 11 ' 5 /i _1 M . Above this mass, / inf w 0.03 
works well across over two decades in Mj n f, consistent with 
results from spatial clustering. 

Note that the satellite fraction remains non-zero up to 
Mi n f ~ 10 14 fe~ 1 M ( y), resulting from low redshift cluster- 
cluster mergers (e.g. lCohn fc Whitdl2005h . This implies that 
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~ 10% of galaxy clusters should be expected to host a 
satellite BCG, as supported by observations whic h shown a 
numb er of clusters hosting BCG-BCG pairs (e.g.. iLiu et al.l 
l2009h . In a sample of rich clus ters from the SPSS MaxBCG 
catalog, 15% host two BCGs ijKoester et al.ll2007h . 

We can also use our analytic model for satellite removal 
]to examine how different Cd yn values for dynam- 



ical friction affect the satellite fracti on, as shown in Fig . [8] 



(bottom). Discounting the results of iTinker et al.l (|2007f ) at 
low mass, which exhibit a flattening that is not seen in any of 
the other samples, our simple analytic model provides good 
agreement with observations using Cd yn = 0.2 — 0.3. This 
provides a consistent picture, since using Cd yn = 0.2 — 0.3 
matches the subhalo HOD using / inf = 0.01 - 0.03 (Fig. 0), 
which is also the range of /i n f that best matches observed 
clustering. Additionally, this analytic model works well even 
below our mass threshold for reliably tracking satellite sub- 
halos, showing no rollover at low mass since it depends only 
on resolving halos prior to infall. 

The dotted curve in Fig. [H] (bottom) shows the satellite 
fraction if satellites are never removed (merge with central 
or tidally disrupt), which we measure by analytically re- 
taining every infalling halo, as in Fig. [5] The slight deficit 
compared with the full subhalo catalog at high mass arises 
because satellites in this analytic model do not gain mass 
via merging with other satellites. While this "no removal" 
scenario provides a reasonable estimate at high mass, where 
halos have formed only recently and so their average times 
since infall are short, it grossly overpredicts the satellite frac- 
tion at lower mass. This demonstrates the clear failure of a 
scenario in which satellite galaxies never merge/disrupt. 

Overall, we find that Cdyn ~ 0.25 both matches well to 
the satellite removal criteria that agrees with spatial cluster- 
ing, and it directly produces sate l lite fr actions in agreement 
with observations. IConrov et al. (2007) apply a similar dy- 
namical infall model as ours to halo merger trees to examine 
the growth of the Luminous Red Galaxy (LRG) population. 
They find that Cd yn = 0.1 fit observations well, though they 
were unable to put precise constraints on this value. Fig- 
ure [5] (bottom) shows that such a value would underpredict 
the observed satellite fraction across a broad mass range. 

Cd V n = 0.2 also agrees well with a similar analysis by 
(I Wetzel et al.ll20 09a). who match the evolution of satellite 
subhalos in simulations at z > 1. They also pointed out 
that, given typical satellite orbital circu larity distributions 
ijZentner et al.l 120051 : Ijiang et al.l l2008h. the detailed fits 
to dy namic al friction in f all tim es of iBovlan-Kolchin et al.l 
(120081 1 and Ijiang et alj (|200Sl) predict C dyn = 0.06 and 
0.14, respectively. The latter value, obtained from an anal- 
ysis of a cosmological simulation not dissimilar from ours, 
is marginally consistent with our results. However, the for- 
mer clearly is not. This can partially be attributed to 
the different mass d e pende nce in the parametrization of 
Boyl an-Kolchin et al.l 1)20081) : they used a fit similar to Eq.[2] 
but with (Mhaio/Msat.inf ) in the numerator. Their added 
exponent factor doubles their merging/disruption time (with 
respect to ours) for a 10 : 1 mass ratio merger. While this 
does provide better agreement, it still leads to a 10 : 1 in- 
fall timescale that is half of ours. Given that their fit was 
obtained from an analysis of significantly higher resolution 
mergers of two isolated, spherical (Hernquist profile) halos, 



this highlights the importance of realistic cosmological set- 
tings for calibrating satellite merging/disruption timescales. 

Finally, we note that the above analytic model of satel- 
lite removal is dependent on halo definition. In particular, if 
halos are defined to have a larger radius, then satellites will 
accrete sooner and hence the analytic removal time will be 
longer. As compared with the commonly used FoF(& = 0.2), 



Tb=0. 168 



0.86r{, = o.2. Thus, the best-match Cd y 



0.25 



values represents a lower limit to that of satellites in halos 
using FoF(fe = 0.2)0 



4.3 Cluster Satellite Luminosity Function 

Since the correlation function and satellite fraction conflate 
satellites within a variety of host halo masses, we finally test 
our subhalo catalog by examining satellites specifically in the 
densest environments: galaxy clusters. We compare the lu- 
minosity function of satellites within galaxy clusters in our 
simulati on with that determ ined from the SDSS MaxBCG 
catalog llHansen et al . 2009). Here, we abundance match our 
subhalo catalog at z = 0.25 to the SDSS i-band galaxy lumi- 



ixy J 

nosity function at that reds hift (ISheldon et al J 120091) . from 



which the cluster catalog in IHansen et alT l2009) is drawn 



Additionally, we assume 0.15 dex scatter in luminosity at 
fixed subhalo mass. We then compute the luminosity func- 
tion of satellites i n clusters of a given r ichness, where we de- 
fine richness as in IHansen et al.l (|2009f ) as the number of red 
sequence satellite galaxies with L > 0.4L* within the clus- 
ter virial radius. We use their fit to the red galaxy fraction 
vs. cluster richness (their Eq. 13) to scale our total satellite 
counts to those of red galaxies (this fraction is 65 — 80% in 
the richness range we probe). 

Figure [5] shows the results of our subhalo catalog, us- 
ing a remova l thres hold of / inf > 0.02, as compared with 
IHansen et al. ( 2009). We find good agreement for both high 
and low bins of cluster richness. Since Mi = — 19 corresponds 
to a subhalo infall mass of 10 11 ' 4 h~ 1 M ( 7 ) , the agreement in 
Fig. [Sjacross a wide range of luminosities demonstrates that 
our subhalo catalog, with a model for merging/disruption, 
successfully traces the galaxy population in the densest envi- 
ronments down to our expected numerical resolution thresh- 
old. Figure [5] does exhibit some excess of very bright cluster 
satellites in our catalog, though the limited number of mas- 
sive clusters in our simulation volume precludes a definitive 
test in that regime. 



5 HIGH REDSHIFT 

5.1 HOD and Analytic Removal 

Examining the effects of removal thresholds on the HOD at 
higher redshift, we find little change in the fractional reduc- 
tion in the satellite HOD as /i n f is increased. This supports a 



4 In principle, changing halo definition can also change the mass 
ratio in Eq. [2] and hence the removal time, since halo concentra- 
tion scales weakly with mass. However, the mass enhancement 
from FoF(fe = 0.168) to FoF(6 = 0.2) is nearly independent of 
concentration and redshift. The relation to spherica l overdensity 
halos will depend on concentration and redshift (see lWhitell200ll , 
for more on halo conversion). 
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Figure 9. Luminosity function of satellites in galaxy clusters 
at z = 0.25, for lower (top) and higher (bottom) cluster richness 
bins. Subhalos a re abundance match ed to the SDSS i-band lumi- 
nosity function (Shel don et alJfeoOSl) . assuming 0.15 dex scatter 
in luminosity at fixed subhalo mass, and satellites retained if their 
/inf > 0-02. S o lid cu rve shows fit to SDSS MaxBCG catalog from 
lHansen et ail ||2009|) . 



picture where satellite subhalo mass stripping is quite rapid 
after infall, such that the shorter average times since infall 
at higher redshift do not lead to a significantly less stripped 
satellite subhalo population. 

We also examine our analytic model for satellite removal 
timescale (Eq.0) at higher redshift. In particular, we test the 
validity of the scaling of td yn with f Hubble- Figure [TOl shows 
the HOD at z = 1 and 2 for both the subhalo catalog and 
analytic model, using the same /i n f and Cdyn values as in 
Fig. [5] (bottom). The agreement at all host halo and satellite 
masses remains relatively robust given the simplicity of the 
model. This agreement is partially a result of the decreas- 
ing dependence of the analytical model on Cd yn at higher z, 
when halos of a given mass have formed and been accreted 
more recently, so allowing longer infall times does not sig- 
nificantly increase the number of satellites. This is demon- 
strated by decreasing difference in the HOD for Cdyn = 0.2 
(dashed curves) as opposed to the "no removal" scenario 
(dotted curves). Moreover, using Cdyn > 0.3 produces no 
additional enhancement in the HOD at z > 2. 
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Figure 10. Same as Fig. [5] but at z = 1 (top) and 2 = 2 (bot- 
tom). A given Cd yn value maps nearly as well to a given removal 
threshold, /j n f, at high z as at z = 0.1. Dotted curves show keep- 
ing all satellites ever accreted. At higher z, this case exhibits less 
excess with respect to the dynamical infall prescription. 



This result has clear implications for modelling galaxy 
evolution at high redshift. Since the details in modeling 
satellite galaxy infall matter less at higher redshift, nearly 
all uncertainty in understanding galaxy evolution at high 
redshift lies in the gas/radiation physics and not in the dy- 
namics of dark matter. 



5.2 Spatial Clustering at z ~ 1 

We also examine whether our criteria for satellite removal 
hold at higher redshift by comparing with galaxy clustering 
results at z ~ 1 from the De ep Extragalactic Evolution- 
ary Probe 2 (DEEP2) survey (jCoil et alj|2006h . as shown 
in Fig. 1111 Here, we use abundance matching to match 
the number densities of B-band luminosity threshold galaxy 
samples in DEEP2. We find that using abundance matching 
with no scatter in the Lb — M\ n i relation causes our subhalo 
catalog to exhibit significant excess in clustering as com- 
pared with DEEP2 for all but the lowest few r p bins. We 
find improved agreement using significant (0.6 dex) scat- 
ter, as shown for two luminosity thresholds corresponding 
to subhalos above our robust resolution limit. Using even 



© 2009 RAS, MNRAS 000.ITHT91 



What determines satellite galaxy disruption? 13 




i i i i i 1 1 1 1 1 i i i i i i 1 1 1 

0.1 1 10 

r p [h-'Mpc] 



Figure 11. Projected auto-correlation function for several re- 
moval thresh olds, /j n f, as com pared with the observed clustering 
at z ~ 1 of ICoil et a"D (|2006h . M B < -19.5, -20.5 samples cor- 
respond to subhalo M inf > 10 115 , 10 12 /i _1 M . Abundance 
matching here assumes 0.6 dex scatter in Lb at fixed M; n f. 



larger scatter improves the agreement on large scales, but it 
underpredicts clustering on smaller scales. 

At z ~ 1, the small scale (r p < lh^ 1 Mpc) clustering 
provides a consistent picture with the results at z = 0.1, 
favoring a removal threshold of /i n f = 0.01 — 0.05. However, 
despite the large luminosity-mass scatter, our subhalo cata- 
log overpredicts the clustering as compared with DEEP2 on 
large scales. Curiously, at face value, DEEP2 clustering at 
r p ~ 10 ft -1 Mpc favors a model with no satellite galaxies. 

There are a number of possible reasons for why 
DEEP2 suggests increased scatter and/or reduced cluster- 
ing. DEEP2 selects on B-band luminosity, which is more sus- 
ceptible to recent star formation, and since satel lites at z ~ 1 
are redder than centrals at a given stellar mass (|Zheng et al.l 
120071 : iTinker fe Wetzelll2009h . this could bias the clustering 
low since satellites preferentially live in more massive halos, 
which are more highly biased. As discussed in £12.21 it is pos- 
sible that the L — Mi n f relation has more scatter at higher 
redshift. However, large scatter alone is unable to match the 
clustering at all scales. Clustering in DEEP2 could be biased 
low because their selection criteria miss ~ 10% of red galax- 
ies in their survey volume, although this effect is likely to be 
mild both because red galaxies form the minority of the pop- 
ulation at all luminosity thresholds and because w p (r p ) is 
nearly the same f or red and blue galaxies at r p ~ 10 h -1 Mpc 
l|Coil et al.ll2008l ). Alternately, if DEEP2 is missing an ap- 
preciable number of galaxies in their survey volume (regard- 
less of color), our abundance matching method will corre- 
late to the observed galaxies artificially high subhalo masses, 
leading to enhanced clustering. Or the disagreement could 
be a statistical fluctuation. The DEEP2 survey volume of 
~(100/i _1 Mpc) 3 is susceptible to sampling variance, and 
the clustering at different scales are highly correlated, so 



the discrepancy may not be statistically significant. Since 
we do not have a covariance matrix for the observations, we 
are unable to det ermine the goodness- of-fit. 

Interestingly, IConrov et al.l (|2006l) find good agreement 
with DEEP2 clustering using similar method based on 
abundance matching against subhalo catalogs. Their better 
agreement (lower clustering) may stem from their different 
assumed cosmology (as = 0.9), and/or their simulation box 
sizes (80 and 120 h~ Mpc), significantly smaller than the 
one used here. Suggestively, a recent HOD analysis using 
color information finds a similar clustering trend as ours, 
namely, good agreemen t at small scales but e xcess cluster- 
ing at r p > 1 h' 1 Mpc l|Tinker fc Wetzelll2009h . 

Thus, while our satellite removal criteria are consistent 
with clustering at z ~ 1 in DEEP2 at r p < Ih^ 1 Mpc 
given Lb — Mi n f scatter, the above uncertainties and lack of 
agreement in large-scale clustering do not allow us to obtain 
strong constraints from these higher redshift data. 

5.3 Evolution of the Satellite Fraction 

Finally, we explore how the satellite fraction evolves with 
redshift. Figure [12] (top) shows the satellite fraction as a 
function of Minf at various redshifts, using /j n f = 0.02. At 
a fixed infall mass, the satellite fraction monotonically in- 
creases with time, arising from the decreasing dynamical 
friction in fall time with respe ct to the typical halo merger 
timescale l|Wetzel et al.l20"09a[ ) . The evolution of the rollover 
at low mass, arising from numerical disruption, indicates 
that the resolution limit decreases with higher redshift, when 
satellites have been accreted more recently. 

Since the M sta r — Mi n f relation evolves with time, Fig. 1121 
(bottom) shows the satellite fraction instead as a function 
of M Btar , obtained b y abundance m atching against the SMF 
of ICole et all (|200ll) at z = 0.1 and iMarchesini et al.l (|2009h 
at higher redshift. Points show the satellite fraction at fixed 
Mi n f values. At high mass, the M atar — Mi n f relation exhibits 
only mild evolution, but at lower mass, subhalos at a fixed 
Minf grow rapidly in M star . 

This result has two important implications. First, in 
examining a galaxy population of fixed stellar mass, simu- 
lation resolution becomes even less of a limitation at higher 
redshift than Fig. [T2] (top) would indicate, since at the low 
mass end galaxies of a fixed Af star move to a higher M star 
with redshift. This is evidenced by the lack of rollover in the 
stellar mass satellite fraction at higher redshift. 

Second, the growth of M 3 t ar at fixed M in f for M star < 
1O 11 M0 demonstrates how SHAM implicitly allows for satel- 
lite star formation. For example, if a satellite were accreted 
at z = 1 with M ata ,r = 4 x 1O 1O M0, following the points of 
fixed Minf in Fig. [12] (bottom) shows that SHAM-assigned 
stellar mass would increase by ~ 20% by z = 0.1. 



6 COMPARISONS WITH OTHER REMOVAL 
CRITERIA 

As discussed in the introduction, various authors have used 
different criteria to define satellite galaxy removal via dark 
matter subhalos. Here, we examine how other criteria com- 
pare against our fiducial case, the subhalo bound mass to 
infall mass ratio, /inf. 
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Figure 12. Satellite fraction vs. mass at various redshifts, us- 
ing a removal threshold of /; n f = 0.02. Top: Subhalos binned 
in Mj n f. The turnover in satellite fraction occurs at lower 
M; n f at higher z, indicating that resolution effects less crit- 
ical at higher z. Bottom: Subhalos binned in M s t ar , based 
on abundance matching. Points indicate values at fixed M; n f 
of 1O 1175 ,1O 12o ,1O 12 - 5 ,1O 13 - o /i- 1 M , highlight the evolving 
M s tar — M inf relation. 



As shown by our full tracking scheme (keeping subhalos 
to 50 particles), using a fixed mass criterion can reasonably 
match observations if the mass threshold corresponds to 1 — 
5% of the infall mass. As simulations increase in resolution, 
it thus becomes increasingly necessary to model for satellite 
removal before subhalos reach numerical disruption. 

Models which remove a satellite subhalo after its mass 
has been stripped to a value less than the m ass within a fixed 
fract i on (adis) of its NFW scale radius (e.g., Taylor fc Babull 
120041 : IZentner et"ai1l2005l : IMaccio et al.ll2009T l are largely de- 
generate with our fiducial /mf criterion, since subhalo con- 
centration scales only weakly with mass. For a typical halo 
concentration of c = 10, our results suggest adis ~ 0.3, 
within the range 0.1 — 1 often assumed. 

Given th at total baryon i c to d ark mass ratio in a halo 
is nearly 10%. IStewart et all |2009h assumed that a satellite 
merges with its central galaxy when it drops below 10% 
of infall mass. However, we find that /i n f = 0.1 noticeably 
underpredicts the satellite population. A lower threshold is 
needed, likely both because the stellar to dark mass ratio is 



significantly lower (~ 4%), and because any baryonic mass 
will be more compact than the 10% most bound dark mass. 

Since the M s tar/Afj n f ratio varies with subhalo mass 
(Fig.[T]), it is not obvious that our best-fit removal threshold, 
based on subhalos of M inf > 10 11 ' 5 ft _1 M s (L > 0.2L„), 
can be simply extrapolated to modelling dwarf-spheroidal 
systems around the Milky Way. However, our results are 
consistent with allowed dark mas s stripping fractions to 
match the Milky W ay population (|Penarrubia et all 120081 : 
IMaccio et ai1l2009ri . 

Instead of mass stripping-based criteria, we also ex- 
amine criteria based on satellite tidal radius or distance 
from halo center/central galaxy. To obtain precise central- 
satellite separations, we linearly extrapolate satellite po- 
sitions between subsequent outputs to find the minimum 
separation. For succinct comparison, we compare satellite 
fractions against /j n f = 0.02, the fiducial criterion which 
matches well to observations as shown above. 

As was shown in Fig. [4] a correlation exists between a 
satellite's distance from its central galaxy and falling below a 
given /j n f value. Thus, Fig.[T3](top) shows satellites removed 
when coming within a fixed physical distance from the cen- 
tral galaxy. Despite the stripping-distance correlation, com- 
pared with our fiducial case this model underpredicts satel- 
lites at low mass and overpredicts at high mass. Combined 
with the results from our full tracking model (fixed mass 
threshold), this implies that any successful removal criteria 
must scale with the mass or radius of the satellite subhalo. 

In this vein, the middle panel shows satellite removal if 
a satellite subhalo's tidal (Jacobi) radius falls below a given 
value, where 



r t id 



.4 



1 1/3 



3M halo (< r) 



(4) 



and where r is the physical separation from satellite to 
halo center, M sat is the instantaneous bound satellite mass, 
Mh a io(< r) is the halo mass within r, and A is a con- 
stant that accounts for t he host halo's density profile (see 
iBinnev fc Tremaindl2008l . Eq. 8.107, 8.108). Assuming that 
a satellite disrupts when within its host halo's inner NFW 
profile lets us set A — 1.4, and we measure M sa t and M na i (< 
r) directly from our simulation. As a simple analytic alter- 
native, we have investigated mod elling the halos at infall as 
isothermal spheres (as applied to lBinnev fc Tremainei[2008l . 
Eq. 9.86), which produces similar results. Figure [T3] (mid- 
dle) shows the satellite fraction for tidal radii motivated by 
typical galaxy sizes. Similar to using fixed satellite distance, 
this method predicts too shallow of a dependence on satellite 
infall mass, likely a manifestation of its weak (1/3 power) 
scaling with mass ratio. 

Finally, we consider a removal criterion defined simply 
by a linear scaling relation between satellite halo - host halo 
mass and radius at infall. Specifically, removal occurs when 



D i MTsat.inf 1 n 

r < B -— itha.lo.vir 

V M ha l, 



(5) 



where B is some constant. Note that if satellite and host ha- 
los are isothermal spheres, such that M(< r)/M v i r — r/i? v i r , 
and B = 1, the above criterion is met when the halo mass 
interior to r equals the satellite's infall mass. Figure[T3] (bot- 
tom) shows the results using B — 1 and 0.5. This model 
provides good agreement with our fiducial case because of 
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Figure 13. Satellite fraction vs. subhalo infall mass, using vari- 
ous removal criteria. Top: Satellite removed when it passes within 
a given physical distance of its central galaxy. Middle: Satellite 
removed when its tidal radius falls below the given value. Bot- 
tom: Satellite removed when it passes within radius at which both 
satellite subhalo and host halo have same interior mass (ignores 
satellite mass stripping), both using NFW profiles from simula- 
tion and approximating halos as isothermal spheres. 



the strong (linear) scaling with mass ratio, with B = 0.5 
well matching the overall amplitude. Thus, this method may 
provide a useful approximation for simple analytic models 
of satellite evolution. 

Various works using semi-analytic models grafted onto 
subhalo populations use an analytic model to track satel- 
lite galaxies after subhalo disruption in o rder to match 
small - scale galaxy clustering at z ~ (e.g.. ISpringel et all 
l200ll ; iKitzbichler fc White! I200I ISaro et al.l l2008h . While 
this method does provide a means to overcome limited mass 
resolution, its fundamental assumption, that all satellite 
galaxies eventually merge with the central galaxy, is not 
consistent with the dynamics of subhalos at removal that 
we find, or with observ ed ICL fractions . 

As an alternative, iHenriaues et all (|2008h . using semi- 
analytic catalogs from the Millennium simulation, consider 
all satellite galaxies at 2 = whose subhalo has fallen below 
numerical resolution as tidally disrupting into ICL (imply- 
ing that central-satellite mergers do not occur). While this 
represents a highly simplified model, entirely contingent on 
a fixed numerical resolution mass, it does provide better 
agreement with the observed faint end of the galaxy lumi- 
nosity function an d allows for contributio n to ICL consistent 
with observations. ISeek Kim et al.l l)2009l ) use the same semi- 
analytic models and find better agreement with observed 
spatial clustering by assuming that satellite luminosity is 
instantly reduced at infall to roughly the ratio of satellite 
halo to host halo mass. This model implies nearly total stel- 
lar mass stripping of low mass satellites immediately upon 
infall, in disagreement with our method which tracks mass 
loss over time acr oss all satellite masse s self-similarly. 

Interestingly, iMoster et all (|2009i) use SHAM applied 
to a simulation of similar mass and force resolution to ours, 
and they find that they need to analytically track satellite 
galaxies after subhalo disruption to match galaxy cluster- 
ing, in opposition to our findings. This difference may arise 
because their simulation box size is 100 h' 1 Mpc, and as we 
discuss in the Appendix, simulation volumes smaller than 
200 h~ x Mpc tend to underpredict subhalo spatial cluster- 
ing on all scales. 



7 SUMMARY AND CONCLUSION 

We use a high-resolution dark matter-only A-body simula- 
tion of cosmological volume to track halos and their sub- 
structure, examining the fates of satellite galaxies. Under 
the assumption that galaxies reside at the centers of subha- 
los and that stellar mass is related to infall mass, we assign 
stellar masses to subhalos through abundance matching, 
such that we recover the observed stellar mass function. Our 
subhalo catalog intrinsically incorporates satellite-satellite 
mergers and satellite orbiting beyond their halo radius. 
Thus, we examine different criteria for satellite galaxy re- 
moval (tidal disruption or merging with the central galaxy), 
using only information from the dark matter density field, 
and focusing primarily on the criterion based on the sub- 
halo bound mass to infall mass ratio, /i nt = A/bound /Mint. 
We highlight our main results as follows: 

• Raising the threshold for removal, /i nt , causes a reduc- 
tion in the satellite HOD at all masses, such that more halos 
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host only one galaxy. Additionally, higher /i n f leads to a shal- 
lower HOD slope, implying that modelling satellite removal 
is more important in higher mass halos. 

• Disrupted/merged satellites reside primarily in the in- 
ner halo regions, and raising the removal threshold leads to 
a less concentrated satellite density profile. A small fraction 
of satellites are disrupted out to the halo virial radius. 

• Most satellites fall below the removal threshold at a 
halo radius of 40 — 100 kpc, the distance increasing with 
/inf. The orbits of these satellites are only mildly radial and 
nearly half are directed outward from halo center, indicating 
that a substantial fraction of satellites are disrupted into the 
ICL as opposed to merging with the central galaxy. 

• Our subhalo catalog matched in number density to ob- 
served galaxy samples at z = 0.1 reproduces well the pro- 
jected correlation function at all scales and over a large range 
in subhalo masses using /j n f = 0.01 — 0.03. Using reasonable 
(0.2 dex) luminosity-mass scatter does not change our re- 
sults appreciably. 

• Our subhalo catalog also agrees well with observed 
satellite fractions and cluster satellite luminosity functions 
for similar values of /i n f . A scenario in which satellite galax- 
ies never merge or disrupt is not feasible. 

• Methods for satellite removal based on physical sepa- 
ration from central galaxy or tidal radii predict too shallow 
of a dependence on subhalo mass. Better agreement arises 
from methods that scale linearly with subhalo mass/size. 

• A simple analytical model for satellite removal 
timescale, based only on the satellite-halo mass ratio at in- 
fall and scaling with the Hubble time, reproduces well both 
our subhalo catalog up to high redshift, and the observed 
satellite fraction at z — 0.1. The best-fit Cdyn ~ 0.25 im- 
plies a substantially longer removal time than that derived 
from other dynamical friction models. 

We emphasize that our results are based on examin- 
ing subhalo merging/disruption in dark matter-only simula- 
tions, and the best-fit criteria to observations may change 
somewhat for simulations incorporating hydrodynamics and 
accurate galaxy stellar profiles. For example, gas dynamics 
can affect orbits of the subhalos them selves, though this ap- 
pears to be small (| Jiang et al.l 120091) . More importantly, a 
realistic subhalo stellar component is likely to be more com- 
pact and immune to stripping than the dark matter. Recent 
works find some enhancement of satellite subhalo surviv- 
ability in SPH simulations with star formation than in dark 
matter-onl y simulations, though the strength of t his effect 
is unclear (jDolag et al]|2009l ; IWeinberg et alj|2008h . 

The largest source of uncertainty in our method is the 
assignment of stellar mass/luminosity to subhalos. While 
abundance matching to subhalo infall mass does allow for 
"passive" satellite star formation as the M star — Mjnf re- 
lation evolves, it may underpredict the star formation in 
satellite galaxies, particularly at higher redshift. We have in- 
vestigated how additional satellite star formation after infall 
influences our results by increasing all satellite infall masses 
by a fixed fraction. Because our removal criterion is based 
on mass stripping, this is largely degenerate with increasing 
removal threshold. Thus, to the extent that satellites expe- 
rience significant star formation after infall (beyond what is 
implicitly in SHAM), our best-fit removal threshold may rep- 
resent a slightly underestimate. A more detailed prescription 



for star formation quenching can be related to environmen- 
tal dependence of galaxy color, though we save incorporating 
properties beyond stellar mass or luminosity to future work. 

Despite these uncertainties, our best fit mass strip- 
ping criterion for removal provides a consistent picture with 
galaxy stellar mass to dark mass ratios. Abundance match- 
ing at z = 0.1 gives M star /M su b,mf ~ 4% for M- m i ~ 
10 12 ft _1 M (Fig.HJ. This implies that all of the stellar mass 
is only slightly more compact than the most bound dark 
matter of the same mass. Thus, it is unlikely that galaxies 
within subhalos that fall below / n f ~ 0.01 can remain in- 
tact. If satellite galaxies have been stripped of most of their 
gas, then the gravitational dynamics of the stellar material 
will simply follow that of the dark matter, and since the dark 
mass fraction within a galaxy's radius is close to unity, dark 
mass stripping below /j n f ~ 0.01 must correspond to some 
amount of baryonic stripping. While our binary model of 
merging/disruption is clearly an oversimplification, it does 
provides a good match to current observations. 

Finally, our results imply strong constraints on simu- 
lations needed for robust subhalo populations, both in res- 
olution and volume. Taking / n f = 0.01 as a conservative 
limit of our removal threshold and ~ 30 particles as mini- 
mum resolution for tracking subhalos in dense environments, 
this implies that at least ~ 3000 particles are needed at in- 
fall to robustly track satellites, which is consistent with our 
convergence tests on subhalo mass functions at z = 0.1. 
This resolution requirement is relaxed somewhat at high 
redshift (z > 1), where average satellite times since infall 
are shorter so satellite have experienced less stripping. Ad- 
ditionally, our tests on simulation volumes indicate that box 
sizes < 200 h~ x Mpc have considerable difficulty reproduc- 
ing representative halo mass functions and spatial cluster- 
ing. Taken together, these resolution and volume effects ne- 
cessitate a wide dynamical range, using several billion par- 
ticles, for simulations used to investigate galaxy evolution 
and galaxy large-scale structure. 
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APPENDIX A: IMPACT OF SIMULATION SIZE 
AND COSMOLOGY 

Given finite computational capacity, there is always a trade- 
off between simulation resolution and volume. While our 
high-resolution 200 hT 1 Mpc simulation is able to track halo 
substructure with high fidelity, its limited volume poses 
problems for accurately recovering the high mass end of the 
halo mass function and large-scale spatial clustering. This 
is not merely an issue of limited statistics and sample vari- 
ance: by imposing mean density on a limited volume, we 
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are not probing an entirely representative realization of the 
Universe, and the limited volume truncates power on scales 
larger than the box size. 

Comparing the halo statistics of the 200 ft -1 Mpc and 
720 h' 1 Mpc simulations (which have different resolution 
but the same cosmology and halo finder, as described in 
i)2.1[) . we find that the halo mass function at z = 0.1 in 
the 200 h~ l Mpc simulation is 5% lower than that in the 
720 /i" 1 Mpc simulation from 3 x 10 11 h~ 1 M Q (the reso- 
lution limit of the larger simulation) to 3 x 10 14 h -1 M©, 
with no dependence on mass in this interval. This implies 
that the deficit is driven not by resolution effects, but in- 
stead by truncated large-scale power and sample variance. 
Above 3 x W 14 h~ 1 M Q , the 200 h' 1 Mpc simulation ex- 
hibits a significant deficit in halo density (~ 20%), with 
its most massive halo being 8 x 10 14 7i _1 iWo, as opposed to 
2 x 10 15 /i _1 Mq in the larger simulation. Additionally, we 
find an error of several percent in the halo correlation func- 
tion measured at 10 Mpc in the 200 ft -1 Mpc simulation 
as compared with the larger one, and it is significantly trun- 
cated on scales > 0.1 Lbox- 

To test these finite-volume issues on subhalo cluster- 
ing, we create a second subhalo catalog by mapping our 
high-resolution subhalo catalog onto the halo catalog of the 
larger, less resolved simulation. Based on convergence tests 
from simulations of multiple sizes, the larger simulation ac- 
curately reproduces the spatial clustering on the scales of in- 
terest. We randomly match two halos of the same mass from 
the two simulations, and we take subhalos from the halo in 
the high-resolution simulation and paste them onto the halo 
in the larger simulation. We place the central subhalo at the 
potential minimum of the halo, and we place the satellite 
subhalos using a random overall orientation such that their 
positions with respect to all other subhalos in the host halo 
are retained. Since halos of M > 8 x 10 14 /i _1 Mq do not 
have a counterpart in the high-resolution simulation, we use 
fits to the satellite HOD extrapolated to higher halo mass 
to populate these halos in the larger simulation, assigning 
the satellite radial distributions to follow their halo NFW 
profiles with a random phase. To examine the influence of 
cosmology, we also populate our subhalo catalog into our 
500 h" 1 Mpc simulation which uses n = 1.0 and as = 0.9. 

Figure IA"T1 shows the effects of both simulation size and 
cosmology on the correlation function of subhalos corre- 
sponding to M r < —21.0. Subhalos in the larger simula- 
tion with the same cosmology exhibit somewhat enhanced 
w p (r p ) on the largest and smallest scales. The change on 
moderate to small scales arises primarily because our re- 
population method ignores any environmental effects on the 
satellites, like halo alignment. However, the large scale clus- 
tering amplitude is little affected. We thus conclude that the 
200 h" 1 Mpc simulation is not significantly affected by box 
size on these scales. 

We emphasize, however, that the spatial clustering in 
smaller simulation box sizes is significantly affected on these 
scales. We have performed the same comparisons with ob- 
servations as in Fig. I All with a simulation of the same mass 
resolution and subhalo finder as the 200 ft -1 Mpc simula- 
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r p [h-'Mpc] 

Figure Al. Same as Fig. [6] but using a fixed /j n f = 0.01 for 
subhalos directly from the 200 ft -1 Mpc simulation, and for sub- 
halos populated into host halos in the 720 h~ 1 Mpc simulation 
(same cosmology) and 500 ft - 1 Mpc simulation (different cosmol- 
ogy). w p (r p ) (top) and ratio of simulations to observed w p (r p ) 
(bottom). Also shown is the reduced x 2 °f the fit to observation 
for each simulation. 

tion, but was a smaller 125 h" 1 Mpc boj0 We find a notice- 
able deficit in w p (r p ) at all subhalo masses both on large 
and small scales as compared with observations and our 
200 h^ 1 Mpc box. The former is readily expected, driven 
by mean density and a truncation of large-scale power. The 
latter effect is more subtle, arising from the truncation of 
the halo mass function at high mass, which strongly affects 
the overall satellite population since these halos host many 
satellites. Thus, we stress that in comparing the spatial clus- 
tering of subhalos in simulations to observed galaxies, both 
resolution and finite volume effects need be considered. 

Figure [ATI also shows that changing cosmology to n = 
1.0 and erg = 0.9 leads to a similar but stronger enhancement 
in w p (r p ) which is less consistent with observations. Thus, if 
erg is considerably larger than our fiducial value of 0.8, our 
best-fit /i n f of 0.01 might be an underestimate. 



REFERENCES 

Baugh C. M., 2006, Reports on Progress in Physics, 69, 
3101 

Berlind A. A., Weinberg D. H., 2002, ApJ, 575, 587 
Berrier J. C, Bullock J. S., Barton E. J., Guenther H. D., 

Zentner A. R., Wechsler R. H., 2006, ApJ, 652, 56 
Bildfell C, Hoekstra H., Babul A., Mahdavi A., 2008, MN- 

RAS, 389, 1637 
Binney J., Tremaine S., 2008, Galactic Dynamics: Second 

Edition. Princeton University Press, Princeton, NJ 

5 The initial conditions for this simulation did had a some low 
modes at low k. 



© 2009 RAS, MNRAS 000.011191 



18 Wetzel & White 



Blumenthal G. R., Faber S. M., Flores R., Primack J. R., 

1986, ApJ, 301, 27 
Boylan-Kolchin M., Ma C.-P, 2007, MNRAS, 374, 1227 
Boylan-Kolchin M., Ma C.-R, Quataert E., 2008, MNRAS, 

383, 93 

Brainerd T. G., 2005, ApJ, 628, L101 

Cohn J. D., White M., 2005, Astroparticle Physics, 24, 316 
Coil A. L., Newman J. A., Cooper M. C, Davis M., Faber 

S. M., Koo D. C, Willmer C. N. A., 2006, ApJ, 644, 671 
Coil A. L., Newman J. A., Croton D., Cooper M. C, Davis 

M., Faber S. M., Gerke B. F., Koo D. C, Padmanabhan 

N., Wechsler R. H., Weiner B. J., 2008, ApJ, 672, 153 
Cole S., Norberg P., Baugh et al., 2001, MNRAS, 326, 255 
Conroy C, Ho S., White M., 2007, MNRAS, 379, 1491 
Conroy C, Wechsler R. H., Kravtsov A. V., 2006, ApJ, 

647, 201 

Conroy C, Wechsler R. H., Kravtsov A. V., 2007, ApJ, 
668, 826 

Cooray A., Sheth R., 2002, Phys. Rep., 372, 1 
Davis M., Efstathiou G., Frenk C. S., White S. D. M., 1985, 
ApJ, 292, 371 

Diaferio A., Kauffmann G., Balogh M. L., White S. D. M., 

Schade D., Ellingson E., 2001, MNRAS, 323, 999 
Diemand J., Kuhlen M., Madau P., 2007, ApJ, 667, 859 
Diemand J., Moore B., Stadel J., 2004, MNRAS, 352, 535 
Dolag K., Borgani S., Murante G., Springel V., 2009, MN- 
RAS, 399, 497 
Dubinski J., 1994, ApJ, 431, 617 

Dunkley J., Komatsu E., Nolta M. R., Spergel D. N., Lar- 
son D., Hinshaw G., Page L., Bennett C. L., Gold B., 
Jarosik N., Weiland J. L., Halpern M., Hill R. S., Kogut 
A., Limon M., Meyer S. S., Tucker G. S., Wollack E., 
Wright E. L., 2009, ApJS, 180, 306 

Faltenbacher A., Jing Y. P., Li C, Mao S., Mo H. J., 
Pasquali A., van den Bosch F. C, 2008, ApJ, 675, 146 

Gao L., De Lucia G., White S. D. M., Jenkins A., 2004, 
MNRAS, 352, LI 

Ghigna S., Moore B., Governato F., Lake G., Quinn T., 
Stadel J., 1998, MNRAS, 300, 146 

Gonzalez A. H., Zaritsky D., Zabludoff A. I., 2007, ApJ, 
666, 147 

Gott J. R. I., Juric M., Schlegel D., Hoyle F., Vogeley M., 
Tegmark M., Bahcall N., Brinkmann J., 2005, ApJ, 624, 
463 

Hansen S. M., McKay T. A., Wechsler R. H., Annis J., 
Sheldon E. S., Kimball A., 2005, ApJ, 633, 122 

Hansen S. M., Sheldon E. S., Wechsler R. H., Koester B. P., 
2009, ApJ, 699, 1333 

Hausman M. A., Ostriker J. P., 1978, ApJ, 224, 320 

Henriques B. M., Bertone S., Thomas P. A., 2008, MNRAS, 
383, 1649 

Jiang C. Y., Jing Y. P., Faltenbacher A., Lin W. P., Li C, 

2008, ApJ, 675, 1095 

Jiang C. Y., Jing Y. P., Lin W. P., 2009, ArXiv e-prints 
Kang X., van den Bosch F. C, 2008, ApJ, 676, L101 
Keres D., Katz N., Fardal M., Dave R., Weinberg D. H., 

2009, MNRAS, 395, 160 

Kitzbichler M. G., White S. D. M., 2008, MNRAS, 391, 
1489 

Klypin A., Gottlober S., Kravtsov A. V., Khokhlov A. M., 
1999, ApJ, 516, 530 



Koester B. P., McKay T. A., Annis J., et al., 2007, ApJ, 
660, 239 

Koester B. P., McKay T. A., Annis J., Wechsler R. H., 

Evrard A. E., Rozo E., Bleem L., Sheldon E. S., Johnston 

D., 2007, ApJ, 660, 221 
Kravtsov A. V., Berlind A. A., Wechsler R. H., Klypin 

A. A., Gottlober S., Allgood B., Primack J. R., 2004, ApJ, 

609, 35 

Kravtsov A. V., Gnedin O. Y., Klypin A. A., 2004, ApJ, 
609, 482 

Krick J. E., Bernstein R. A., 2007, AJ, 134, 466 
Kroupa P., Theis C, Boily C. M., 2005, A&A, 431, 517 
Libeskind N. I., Frenk C. S., Cole S., Jenkins A., Helly 

J. C, 2009, MNRAS, 399, 550 
Limousin M., Kneib J. P., Bardeau S., Natarajan P., Czoske 

O., Smail I., Ebeling H., Smith G. P., 2007, A&A, 461, 881 
Limousin M., Sommer-Larsen J., Natarajan P., Milvang- 

Jensen B., 2009, ApJ, 696, 1771 
Lin Y.-T., Mohr J. J., 2004, ApJ, 617, 879 
Lin Y.-T., Mohr J. J., Stanford S. A., 2004, ApJ, 610, 745 
Liu F. S., Mao S., Deng Z. G., Xia X. Y., Wen Z. L., 2009, 

MNRAS, 396, 2003 
Maccio A. V., Kang X., Fontanot F., Somerville R. S., Ko- 

posov S., Monaco P., 2009, MNRAS, pp 1907-+ 
Maciejewski M., Colombi S., Springel V., Alard C, Bouchet 

F. R., 2009, MNRAS, 396, 1329 
Mandelbaum R., Seljak U., Kauffmann G., Hirata C. M., 

Brinkmann J., 2006, MNRAS, 368, 715 
Marchesini D., van Dokkum P. G., Forster Schreiber N. M., 

Franx M., Labbe I., Wuyts S., 2009, ApJ, 701, 1765 
Metz M., Kroupa P., Jerjen H., 2009, MNRAS, 394, 2223 
Mihos J. C, Harding P., Feldmeier J., Morrison H., 2005, 

ApJ, 631, L41 

Mo H. J., Mao S., White S. D. M., 1998, MNRAS, 295, 319 
Monaco P., Murante G., Borgani S., Fontanot F., 2006, 
ApJ, 652, L89 

Moore B., Ghigna S., Governato F., Lake G., Quinn T., 

Stadel J., Tozzi P., 1999, ApJ, 524, L19 
More S., van den Bosch F. C, Cacciato M., Mo H. J., Yang 

X., Li R., 2009, MNRAS, 392, 801 
Moster B. P., Somerville R. S., Maulbetsch C, van den 

Bosch F. C, Maccio' A. V., Naab T., Oser L., 2009, ArXiv 

e-prints 

Murante G., Arnaboldi M., Gerhard O., Borgani S., Cheng 
L. M., Diaferio A., Dolag K., Moscardini L., Tormen G., 
Tornatore L., Tozzi P., 2004, ApJ, 607, L83 

Muzzin A., Yee H. K. C, Hall P. B., Ellingson E., Lin H., 
2007, ApJ, 659, 1106 

Nagai D., Kravtsov A. V., 2005, ApJ, 618, 557 

Natarajan P., De Lucia G., Springel V., 2007, MNRAS, 
376, 180 

Natarajan P., Kneib J., Smail I., Treu T., Ellis R., Moran 
S., Limousin M., Czoske O., 2009, ApJ, 693, 970 

Navarro J. F., Frenk C. S., White S. D. M., 1996, ApJ, 462, 
563 

Ostriker J. P., Tremaine S. D., 1975, ApJ, 202, L113 

Padmanabhan N., White M., Eisenstein D. J., 2007, MN- 
RAS, 376, 1702 

Penarrubia J., Navarro J. F., McConnachie A. W., 2008, 
ApJ, 673, 226 

Peacock J. A., Smith R. E., 2000, MNRAS, 318, 1144 



© 2009 RAS, MNRAS OOO.mff9l 



What determines satellite galaxy disruption? 19 



Purcell C. W., Bullock J. S., Zentner A. R., 2007, ApJ, 
666, 20 

Reichardt C. L., Ade P. A. R., Bock J. J., Bond J. R., 
Brevik J. A., Contaldi C. R., Daub M. D., Dempsey J. T., 
Goldstein J. H., Holzapfel W. L., Kuo C. L., Lange A. E., 
Lueker M., Newcomb M., Peterson J. B., Ruhl J., Runyan 
M. C, Staniszewski Z., 2009, ApJ, 694, 1200 
Rudick C. S., Mihos J. C, McBride C, 2006, ApJ, 648, 936 
Saro A., De Lucia G., Dolag K., Borgani S., 2008, MNRAS, 
391, 565 

Seek Kim H., Baugh C. M., Cole S., Frenk C. S., Benson 

A. J., 2009, ArXiv e-prints 
Seljak U., 2000, MNRAS, 318, 203 

Shankar F., Lapi A., Salucci P., De Zotti G., Danese L., 

2006, ApJ, 643, 14 

Sheldon E. S., Johnston D. E., Masjedi M., Mc Kay T. A., 
Blanton M. R., Scranton R., Wechsler R. H., Koester 

B. P., Hansen S. M., Frieman J. A., Annis J., 2009, ApJ, 
703, 2232 

Simha V., Weinberg D. H., Dave R., Gnedin O. Y., Katz 
N., Keres D., 2009, MNRAS, 399, 650 

Skibba R. A., 2009, MNRAS, 392, 1467 

Smoot G. F., Bennett C. L., Kogut A., Wright E. L., et al., 
1992, ApJ, 396, LI 

Sommer-Larsen J., Romeo A. D., Portinari L., 2005, MN- 
RAS, 357, 478 

Springel V., White S. D. M., Tormen G., Kauffmann G., 

2001, MNRAS, 328, 726 
Stewart K. R., Bullock J. S., Barton E. J., Wechsler R. H., 

2009, ApJ, 702, 1005 
Tasitsiomi A., Kravtsov A. V., Wechsler R. H., Primack 

J. R., 2004, ApJ, 614, 533 
Taylor J. E., Babul A., 2004, MNRAS, 348, 811 
Tegmark M., Eisenstein D. J., Strauss M. A., et al., 2006, 

Phys. Rev. D, 74, 123507 
Tinker J. L., Norberg P., Weinberg D. H., Warren M. S., 

2007, ApJ, 659, 877 

Tinker J. L., Wetzel A. R., 2009, ArXiv e-prints 
Tormen G., Diaferio A., Syer D., 1998, MNRAS, 299, 728 
Vale A., Ostriker J. P., 2006, MNRAS, 371, 1173 
van den Bosch F. C, Aquino D., Yang X., Mo H. J., 
Pasquali A., Mcintosh D. H., Weinmann S. M., Kang X., 

2008, MNRAS, 387, 79 

van den Bosch F. C, Yang X., Mo H. J., Norberg P., 2005, 
MNRAS, 356, 1233 

van den Bosch F. C, Yang X., Mo H. J., Weinmann S. M., 
Maccio A. V., More S., Cacciato M., Skibba R., Kang X., 
2007, MNRAS, 376, 841 

Wang L., Li C, Kauffmann G., De Lucia G., 2006, MN- 
RAS, 371, 537 

Wang L., Li C, Kauffmann G., De Lucia G., 2007, MN- 
RAS, 377, 1419 

Weinberg D. H., Colombi S., Dave R., Katz N., 2008, ApJ, 
678, 6 

Weinmann S. M., Kauffmann G., van den Bosch F. C, 
Pasquali A., Mcintosh D. H., Mo H., Yang X., Guo Y., 

2009, MNRAS, 394, 1213 

Wetzel A. R., Cohn J. D., White M., 2009a, MNRAS, 395, 
1376 

Wetzel A. R., Cohn J. D., White M., 2009b, MNRAS, 394, 
2182 

White M., 2001, A&A, 367, 27 



White M., 2002, ApJS, 143, 241 

White M., Zheng Z., Brown M. J. I., Dey A., Jannuzi B. T., 

2007, ApJ, 655, L69 
White S. D. M., Rees M. J., 1978, MNRAS, 183, 341 
Willman B., Governato F., Wadsley J., Quinn T., 2004, 

MNRAS, 355, 159 
Yang X., Mo H. J., van den Bosch F. C, 2008, ApJ, 676, 

248 

Yang X., Mo H. J., van den Bosch F. C, 2009, ApJ, 693, 
830 

Yang X., Mo H. J., van den Bosch F. C, Pasquali A., Li 

C, Barden M., 2007, ApJ, 671, 153 
Yang X., Mo H. J., van den Bosch F. C, Weinmann S. M., 

Li C, Jing Y. P., 2005, MNRAS, 362, 711 
Zehavi I., Zheng Z., Weinberg D. H., et al., 2005, ApJ, 630, 

1 

Zentner A. R., Berlind A. A., Bullock J. S., Kravtsov A. V., 

Wechsler R. FL, 2005, ApJ, 624, 505 
Zheng Z., Coil A. L., Zehavi I., 2007, ApJ, 667, 760 
Zibetti S., White S. D. M., Schneider D. P., Brinkmann J., 

2005, MNRAS, 358, 949 



© 2009 RAS, MNRAS Q00.fTUT9l 



